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Abstract 

The claim that the light quark mass ratio (m^ — m u )/m s can be extracted from the decay 
width ratio T(j]' —* 7r°7r + 7r _ )/r(r/ / — > 7/7r + 7r~) is critically investigated within a U(3) chiral 
unitary framework. The influence of the recent VES data on the if — > r]7r + TT~ decay is 
also discussed. 
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1 Introduction 



The light quark masses m u , m d , m s are fundamental parameters of Quantum Chromo Dynamics 
and ought to be constrained as accurately as possible. The determination of the light quark 
mass ratios has been the goal of a variety of investigations in low-energy hadron physics, see 
e.g. [1-5]. Of particular interest is the quark mass difference rrid — m u which induces isospin 
breaking in QCD. Moreover, the possibility m u = would provide an explanation for the strong 
CP problem. 

An accurate way of extracting m d — m u is given by the isospin- violating decays 77, rj' — > 
7r°7r + 7r~ and 77,77' — > 3n°. While for most processes isospin-violation of the strong interactions 
is masked by electromagnetic effects, these corrections are expected to be small for the three 
pion decays of the 77 and 77' (Sutherland's theorem) [6] which has been confirmed in an effective 
Lagrangian framework [7]. Neglecting electromagnetic corrections the decay amplitudes are 
directly proportional to m d — m u . 

For this reason, it has been claimed in [8] that the branching ratio r = T(i]' — > 7i 7i + 7i^)/T(rj' - 
777r + 7r~) can be utilized in a very simple manner to extract the light quark mass difference 
Tfid — iTiu- To this aim, it is assumed that 

a) the amplitude ^(77' — > 7r°7r + 7r~) is determined by the corresponding amplitude ^(77' — > 
777r + 7r~) via 

A(r]' -> = e A(r]' -> t/tt+tt") (1) 

with e = (a/3/4) (m d — m u ) / (m s — rh) the 7r°-77 mixing angle and rh = (m d + m u )/2. (Note 
that in [8] the difference m s — rh has been approximated by m s in the denominator of e.) 
Eq. (1) implies that the decay 77' — > 7r°7r + 7r~ proceeds entirely via 77' — > riit + -n~ followed 
by 7r°-77 mixing. 

b) both amplitudes are "essentially constant" over phase space (see the remark in front of 
Eq. (19) of Ref. [8]). 

Based on these two assumptions one arrives at the relation 

_ T(rf -> A+vr-) 3/ m d -m u \ 2 

r " r(77' - 777T+7T-) " lib ^16 V ms ) ' 1 ' 

where the factor 16.8 represents the phase space ratio. Comparison with experimental data — 
for which, so far, only an upper limit exists — would then lead to a prediction for the quark 
mass ratio (m d — m u )/(m s — rh) ~ {m d — m u )/m s . The purpose of the present work is to 
critically examine these two assumptions which lead to the simple relation in Eq. (2). Such an 
investigation is very timely in view of the recent and ongoing experimental activities on 77 and 
77' decays at the WASA facility at COSY [9], MAMI-C [10], KLOE at DA$NE [11] and by the 
VES Collaboration [12,13]. 

An appropriate theoretical framework to investigate low-energy hadronic physics is provided 
by chiral perturbation theory (ChPT) [14], the effective field theory of QCD. In ChPT Green 
functions are expanded perturbatively in powers of Goldstone boson masses and small momenta. 
However, final-state interactions in 77 — > 3ir have been shown to be substantial both in a 
complete one-loop calculation in SU(3) ChPT [15] and using a dispersive framework [16, 17]. It 
is hence important to include final-state interactions in a non-perturbative fashion. 
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In rf decays final-state interactions are expected to be even more important due to larger 
phase space and the presence of nearby resonances. In this investigation, we apply the frame- 
work of U(3) chiral effective field theory in combination with a relativistic coupled-channels 
approach developed in [18, 19] in order to calculate the ratio r. Final-state interactions are 
included by deriving s- and p-wave interaction kernels for meson-meson scattering from the 
chiral effective Lagrangian and iterating them in a Bethe-Salpeter equation. The infinite it- 
eration of the chiral effective potentials generates resonances dynamically. Very good overall 
agreement with currently available data on rj, rf decay widths and spectral shapes has been 
achieved in [18,19]. 

In the next section, we will investigate in our approach if both the assumptions "a" and 
"6" are justified. The inclusion of the recent VES data [13] which provide higher statistics on 
the spectral shape of rf — > ^7r + 7r~ than previous experiments is studied in Sec. 3, while Sec. 4 
contains our conclusions. 



2 Validity of the assumptions 

In the following, we will work with the double quark mass ratio Q 2 

q 2 = m s -m m s + m 
m d -m u m d + m u 

instead of the mixing angle e, since the Kaplan-Manohar reparametrization invariance [20] of 
the chiral effective Lagrangian is respected by Q 2 up to chiral order 0(p A ), whereas e receives 
corrections already at 0(p 2 ). Hence, it is preferable to employ Q in phenomeno logical analyses 
in order to suppress the ambiguity stemming from this reparametrization invariance. 

Following Dashen's theorem which asserts equal electromagnetic corrections for pion and 
kaon masses at leading chiral order [21], Q 2 can be expressed in terms of physical meson masses 
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VDashen - 2 „2 _ ~2 , _,2 _ _,2 _ ^ 1 > ' W 

However, there are various investigations on the size of violations to Dashen's theorem with 
(partially contradictory) results for Q in the range of about 21 . . .24 [22]. The 3ir decays of rj 
and rf provide thus a good opportunity to pin down the value of the double quark mass ratio 
Q 2 [16,23]. 

We will first investigate the validity of assumption "a", i.e. we assume that the decay rf — > 
7r°7r + 7r~ proceeds entirely via rj' — > i]7t + tt~ followed by tx q -7] mixing. This implies for the neutral 
decay A(r)' — > 37r°) = 3e A(r)' — > rj^n ). Employing the amplitudes A(r)' — > 77 27r) from the 
approach advocated in [19] — which are in very good agreement with the data given in [24] - 
one can thus predict the decay amplitudes for A(r)' — > 3n) and calculate both the decay widths 
T(i]' — > 7r°7r + 7r~), Fiji' — > 37T°) and the branching ratios r and T2 = — > 37T )/r(i] f — > rjn ^). 

In [19] a least-squares fit to meson-meson scattering phase shifts and rj, rj' decays has been 
performed. One observes four different classes of fits, i.e. clusters, which describe these data 
equally well, but differ in their predictions for yet unmeasured quantities such as the rj' — > 
7r°7r + 7r~ decay width (for which there exists only a weak upper limit) and the Dalitz slope 
parameters of r]' — > 3n. In fact, as we will see in the next section, inclusion of the recent VES 
data for rj' — > r]n + n^ [13] reduces the number of fit clusters to one, but in the current section we 
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Cluster 1 


Cluster 2 


Cluster 3 


Cluster 4 


Exp. 


rV -> 7r u 7T+7r-) [eV] 


69 ±12 


73±9 


141 ±44 


141 ±26 


< 3800 


r [%} 


0.09 ±0.02 


0.09 ±0.02 


0.17±0.06 


0.17±0.03 


< 4.1 


T{r]' -> 3tt°) [eV] 


116 ± 22 


120 ±16 


217±67 


217 ±40 


315 ±78 


r 2 [%] 


0.26 ±0.05 


0.26 ±0.04 


0.47±0.15 


0.47 ±0.08 


0.74 ±0.12 



Table 1: Decay widths and branching ratios in the chiral unitary approach [19] employing 
assumption "a". 





Cluster 1 


Cluster 2 


Cluster 3 


Cluster 4 


Exp. 


T(rf -> 7T u 7r + 7r-) [eV] 


155 ±7 


155 ±7 


153 ±7 


154 ±5 


< 3800 


r [%} 


0.19 


0.19 


0.19 


0.19 


< 4.1 


F(rf -> 3tt u ) [eV] 


238 ±11 


239 ±10 


237± 11 


239 ±6 


315 ±78 


[%] 


0.52 


0.52 


0.52 


0.52 


0.74 ±0.12 



Table 2: Decay widths and branching ratios in the chiral unitary approach employing assump- 
tions "a" and "6". Since in this case the branching ratios only depend on phase space and Q, 
they do not have an error bar. 

take the data given in [24]. We employ furthermore the upper limit of 1.75 % for the branching 
fraction of rj' — > tt°7t + 7t~ as measured by the VES collaboration [12] which is significantly lower 
than the previous upper limit of 5 % [24] . This tighter bound translates to an upper limit of 
3.8 keV for the partial decay width and reduces the upper limit for r from 10% (as quoted by 
the PDG) to 4.1 %. The pertinent results for the four fit clusters of [19] are well below these new 
upper limits and can be utilized without modification. As already reported in [19], the fit to 
the data does not allow for conclusions on the size of violations to Dashen's theorem since Q is 
treated as an input parameter and variations of Q in the range of 20 ... 24 lead to equally good 
fits within our approach. Hence, we will set Q = 24.1 in our calculations — the value predicted 
by Dashen's theorem. The results for the branching ratios obtained by employing assumption 
"a" and Q = 24.1 are shown in Table 1. The ratios are obtained by explicitly performing the 
integration of the amplitudes over phase space. Obviously, assumption "a" is not justified — at 
least for the neutral decay where, in particular, clusters 1 and 2 are in clear disagreement with 
experiment. 

Next, we employ in addition assumption "6". This is achieved by averaging the 77' — > 77 27r 
amplitudes over phase space which are then employed for rj' — > 3n by means of assumption "a". 
The results are displayed in Table 2. One observes that for clusters 1 and 2 the decay widths 
into 3tt and hence the ratios r, r 2 increase, while the changes for clusters 3 and 4 are rather 
moderate. However, recall that the Dalitz plot parameters of the approach [19] clearly indicate 
that the assumption of a constant amplitude is not justified for rf — > 7r°7r + 7r~, particularly for 
clusters 3 and 4. The partial compensation of the effects of assumption "a" in clusters 1, 2 and 
the moderate changes in clusters 3, 4 are therefore purely accidental. 

We conclude that both assumptions "a" and "6" are not justified. This is further substanti- 
ated by comparison of r and r 2 in Tab. 2 with the respective values from the full chiral unitary 
approach shown in Tab. 3. The values are in clear disagreement and, hence, both assumptions 
are not appropriate — at least within the chiral unitary approach. 

Finally, we would like to investigate the differences which result if assumption "a" is replaced 
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Cluster 1 


Cluster 2 


Cluster 3 


Cluster 4 


Exp. 


rV -> 7r u 7T+7r-) [eV] 


470 ± 200 


520 ±200 


740 ±420 


620 ±180 


< 3800 


r [%] 


0.58 ±0.24 


0.66 ±0.27 


0.92 ±0.52 


0.77±0.21 


< 4.1 


T{r]' -> 3tt°) [eV] 


331 ±24 


326 ±28 


330 ±33 


336 ±21 


315 ±78 


r-2 [%] 


0.73 ±0.06 


0.72 ±0.06 


0.71 ±0.07 


0.73 ±0.05 


0.74 ±0.12 



Table 3: Decay widths and branching ratios in the chiral unitary approach [19]. 





Cluster 1 


Cluster 2 


Cluster 3 


Cluster 4 


Exp. 


T(rf -> 7T u 7r+7r-) [eV] 


2450 ± 1930 


1720 ±1160 


260 ±260 


290 ±290 


< 3800 


r [%] 


2.96 ±2.30 


2.10 ±1.40 


0.34 ±0.34 


0.37±0.37 


< 4.1 


F(rf -> 3tt u ) [eV] 


1080 ±840 


800 ±550 


120 ±120 


120 ±120 


315 ±78 


[%] 


2.34 ±1.79 


1.73 ±1.19 


0.28 ±0.28 


0.28 ±0.28 


0.74±0.12 



Table 4: Decay widths and branching ratios in the chiral unitary approach if isospin-breaking 
takes place solely via 7r -?/ mixing. For the fits of Clusters 3 and 4 this mixing angle can 
actually become zero leading to vanishing decay widths — > 3ir) and branching ratios r, r 2 . 



by the decay mechanism where rf — > 3n occurs due to ir°-r)' mixing followed by a (virtual) 
transition ir° — > 3ir. Employing the relation A(rj' — > 3n) = e' A(ir° — > 37r) with e' being the 
7i°-i]' mixing angle [18] we find the values shown in Table 4. Assuming the rf — > 3n decays to 
proceed via this mechanism introduces a huge uncertainty and leads to different ratios r and r 2 . 
This underlines the observation that the decays rf — > 3n cannot be expected to simply proceed 
either via n°-rj or 7T -?/ mixing. In particular, the isospin-breaking transition due to the quark 
mass difference m d — m u cannot be completely assigned to ir -r] mixing as done in assumption 
"a". Despite its appealing simplicity, the crude estimate given in Eq. (2) is certainly not suited 
to precisely determine the double quark mass ratio Q 2 . In fact, even at leading chiral order 
the rf — > 37r decay amplitude is not entirely due to 7i°-i] mixing. There is also a contribution 
from an isospin-violating ?/37r-vertex from the explicit chiral symmetry breaking part of the 
Lagrangian at second chiral order, see e.g. Ref. [18]. 

On the other hand, employing the chiral unitary approach of [19] does not lead to a conclu- 
sive extraction of Q due to the present experimental situation. From a fit to the data in [24] 
supplemented by Dashen's theorem one obtains the decay width ratio r = (0.35 . . . 1.5)% which 
is larger than the value of 0.18% quoted in [8]. Note also that there is a tendency to even larger 
values of r if Q is lowered, e.g., for Q = 22 we obtain the range r = (0.4 . . . 2.8)%. 

3 Inclusion of the VES data for rj' — ► rj7r + 7r~ 

In this section we study the changes in our results that occur if the recent VES data on the 
spectral shape of 7/ — > i]7i + ii~ [13] are taken into account. Note that the most recent analysis 
of the VES Collaboration [13] has not yet been included in Ref. [24]. The VES data have much 
higher statistics on the Dalitz slope parameters than previous experiments and by including 
them in the fit we obtain the results shown in Tab. 5. Since the amplitudes for rj' — > r]7i + ii^ 
and rj' — > i]7i°7T° are equal in the isospin limit and deviations are thus isospin-breaking and 
small in our approach, we only include the leading Dalitz parameter a of v( — > i]7i°7i° [25] and 
omit the higher ones which are — assuming only small isospin-violating contributions — not quite 
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rf 


— > r]7Y + n 






a 


b 


c 


Theory 


-0.109 


-0.087 


-0.036 


Exp. 


-0.127 ±0.024 


-0.106 ±0.042 


-0.082 ±0.025 




i]' 


— > 1]7l°7l° 






a 


b 


c 


Theory 


-0.123 


-0.104 


-0.041 


Exp. 


-0.116 ±0.026 







Table 5: Results for the Dalitz plot parameters of — > rjinr if the VES data are included in 
the fit. 





Theory 


Exp. 


r(V -> vr u 7r + 7r-) [eV] 


3120 


< 3800 


r [%] 


3.9 


< 4.1 


r(?/ -> 3tt u ) [eV] 


330 


315 ±78 


[%] 


0.73 


0.74 ±0.12 



Table 6: Decay widths and branching ratios if the VES data are taken into account. 

compatible with the new results of the VES experiment for rf — > r]ii + 7i^. Our results are in 
good agreement with the Dalitz plot parameters extracted from the VES experiment. In Tab. 5 
only the best least-squares fit is shown which is sufficient to discuss the qualitative changes of 
the results compared to those of Sec. 2. Note also that we have supplemented our fitting routine 
by a conjugate gradient method [26] and hence the numerical values have improved with respect 
to [19]. 

Our results for the rf — > 37r decay widths and width ratios are displayed in Tab. 6. It is 
important to emphasize that the inclusion of the VES data reduces the number of fit clusters 
to one and we observe indeed one global minimum. There is, however, a strong tendency of 
the fits towards the upper limit T(rj' — > 7r°7r + 7r~) < 3.8 keV and, in fact, slightly improved 
fits with a smaller \ 2 value can be obtained if this upper limit is omitted. (In this case, the 
best overall fit leads to the width T(i]' — > 7r°7r + 7r~) = 5.73 keV.) The result for the width 
ratio, r = 3.9%, has thus increased if the VES data are taken into account. Furthermore, the 
amplitude A(r)' — > 7r°7r + 7r~) fluctuates strongly over phase space with slope parameters which 
can be more than one order of magnitude larger in size than those obtained in [19]. The Dalitz 
plot distribution \A(r)' — > 7r°7r + 7r~)| 2 can therefore not be properly described by a low-order 
polynomial in the usual expansion variables x and y, see [19] for definitions, and it can definitely 
not be assumed to be constant over phase space. 

The reason for both the large decay width F(r)' — > 7r°7r + 7r~) and the strong fluctuations over 
phase space are mainly due to a large contribution from isospin 1 = 1 p-waves in the final-state 
interactions of the decay. While for 1 = 1 p-waves the uncharged two-particle channels are 
C-even and, due to C-invariance, do not couple to C-odd channels related to the p°(770) as 
already pointed out in [19], the coupling of charged channels to the p ± (770) is not forbidden. In 
fact, an important feature of the fits including the VES data compared to those without these 
is the large enhancement of the r/'n^ 1 — > 7r°7r ± coupling which also determines the importance of 
the p ± (770) in this decay. The pertinent Dalitz plot is shown in Fig. 1 and exhibits signatures 
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-1.0 -0.5 0.0 0.5 1.0 

x 

Figure 1: Dalitz plot distribution \A(r]' — > 7r°7r + 7r~)| 2 of the best overall fit including the VES 
data for rj' — > i]7i + ii^ [13]. The distribution is normalized to unity at x = y = (see [19] for 
definitions of x, y). The p-wave contributions to 7r 7r + (7r°7r~) rescattering vanish on the rising 
(falling) dashed line and the invariant energies associated with the p ± (770) are indicated by 
the dotted lines. 

of the p ± (770). Note, however, that these resonances do not appear as bands of increased 
amplitude at fixed two-particle energies (dotted lines in Fig. 1), since the p-wave contributions 
have a kinematical zero in the middle of these bands as indicated in Fig. 1 (dashed lines). Thus 
the amplitude only peaks at the edge of the Dalitz plot. Moreover, due to the symmetry of the 
amplitude under 7r + <-> 7r~ exchange (C-invariance) the p + , p~ peaks interfere constructively 
on the symmetry axis producing a pronounced peak structure at the top of the Dalitz plot, 
where the invariant mass of the 7r + 7r" system is minimal. These features of the Dalitz plot of 
a pseudoscalar meson decaying into three pions have been pointed out long ago in [27]. 

4 Conclusions 

In this work, we have critically investigated the claim of Ref. [8] that the light quark mass ratio 
(rrid — m u ) / (m s — m) can be extracted from the decay width ratio r = T(rf — > 7r°7r + 7r~)/T(r/ — > 
i]7T + 7T~). In order to study this issue we have employed a U(3) chiral unitary framework 
developed in [18,19] which is in very good agreement with the 77, rj' data on widths and spectral 
shapes. Our results clearly indicate that the two underlying assumptions of [8] in order to 
arrive at a relation between r and (m^ — m u )/{m s — m), i.e., that a) the decay rj' — > n n + n^ 
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proceeds entirely via the decay rj' — > ??7r + 7r _ followed by 7i -r] mixing and that b) the decay 
amplitudes are constant over phase space, are not justified at all. The results from the full 
chiral unitary approach are in plain disagreement with these two assumptions. Moreover, the 
present experimental situation which is used as input to fit the parameters of the chiral unitary 
approach does not allow for a precise determination of the double quark mass ratio Q 2 from r. 

Inclusion of the recent VES data on the rj' — > rj-K + -n~ spectral shape reduces the uncertainty 
of the fit results to some extent. In this case, the overall fit to 77, rj' data yields for rj' — > it it + it~ 
a large contribution from the isospin 1=1 p-wave in the final-state interactions which can 
be attributed to a large coupling to the p ± (770) resonances while contributions related to the 
p°(770) are forbidden by C-invariance. More precise data on 77 and rj' decays are needed in order 
to eventually clarify this issue. An improvement of the experimental situation is foreseen in the 
near future due to the upcoming data from WASA at COSY [9], MAMI-C [10] and KLOE at 
DA$NE [11]. 
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